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摘 要 土壤 含水 量 是 影响 半 干 旱 区 农作物 生长 的 重要 因素 。 为 了 准确 测定 土壤 含水 量 的 变化 动态 以 指导 农 
业 高 效用 水 ， 近 年来, 利用 地 球 物理 测量 方法 研究 高 分 辨 率 的 水 流入 渗 ， 已 经 越 来 越 受 到 欢迎 和 重视 。 本 文 以 
陇 中 半 干 旱 区 玉米 田 为 例 ， 通 过 在 土壤 表面 布置 电极 ， 利 用 高 密度 电阻 率 成 像 法 (ERT) 对 降雨 前 后 土壤 二 维 剖 
面 进 行 电阻 率 数据 测量 ， 实 现 对 土壤 二 维 剖 面 电阻 率 值 和 含水 量 监测 ， 解 释 不 同 条 件 下 土壤 含水 量变 化 的 原 
因 ， 建 立 陇 中 半 和 干旱 区 农田 土壤 电阻 率 和 含水 量 之 间 的 相关 关系 。 结 果 表 明 : 降水 入 渗 使 得 二 维 剖 面 土壤 电阻 
率 整体 呈 明 显 降低 趋势 ,， 反 演 得 到 的 电阻 率 图 像 局 部 电阻 值 “高 - 低 - 高 ”的 变化 过 程 ， 与 一 次 降水 过 程 前 后 “ 干 - 湿 - 
干 ”* 的 循环 过 程 一 致 。 土 壤 仿 水量 实测 值 与 估计 值 之 间 有 和 较为 显著 的 线性 关系 (R*=0.651 8, n=96)。 在 0~2.0 m 
深度 范围 内 ， 总 体 估计 偏差 较 小 ,为 0.74%; 土壤 含水 率 的 估计 精度 较 高 ,为 2.64%。0~0.5 m 土 层 (H1) 含 水 量 
监测 探头 分 布 密集 ， 数据 采集 较为 准确 ， 故 Hi 层 估计 精度 略 高 于 0.5~2.0 m É (H5). 相 比 之 前 利用 实测 工具 进 
行 野 外 测量 ,ERT 测量 方法 精度 较 高 。 本 文 提供 了 一 个 高 分 辩 率 的 土壤 结构 二 维 分 布 与 水 分 运 移 过 程 的 图 像 
同时 为 实现 精确 和 高 效 的 农业 用 水 管理 提供 一 种 新 途径 。 
关键 词 ”电阻 率 成 像 ”土壤 含水 量 降水 入 渗 ” 非 饱和 土壤 时 域 反射 ” 陇 中 半 干 时 区 
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Abstract 


placed on the soil surface. Then SWC and electrical resistivity were monitored in two-dimensional soil profile in order to 
explain the variations in SWC under different conditions and to determine the correlation between electrical resistivity of each 


soil horizon and its water content in Longzhong semi-arid region. Soil temperature was also discussed in the ERT data 


Soil water content (SWC) is a key factor influencing crop growth in semi-arid regions. Monitoring the variations in 
SWC is critical for agriculture. In recent years, geophysical survey has been used in the study of high-resolution detection of 
water infiltration and it is hugely popular because of its zero-damage to soil micro-structure and easy measurement. With a 
maize field in Longzhong semi-arid region as the case study, we used high-density electrical resistivity tomography (ERT) to 
measure (before and after the precipitation) two-dimensional electrical resistivity of the column of soil below an electrode 
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interpretation. The results revealed the potential of ERT to improve soil and agronomic studies. Vertical distributions of 
two-dimensional electrical resistivity image inverted from measured data were different. On the whole, natural rainfall 
infiltration reduced the trend in two-dimensional electrical resistivity. The trend of local electrical resistivity of the inverted 
ERT images was one of "high-low-high", quite consistent with the *dry-wet-dry" cycle of the precipitation process. Two- 
dimensional sections of SWC calculated using ERT showed a reliable linear correlation (R? — 0.651 8, n = 96) between the 
estimated and measured SWC in the root-zone horizon, with a slope approximately equal to 1. Within the depth range of 0-2.0 m, 
the precision of the calculated specific SWC quantified by the root mean square error (RMSE) was 2.6496, with a bias 
corresponding to an overestimation of 0.74%. The densely distributed SWC detectors installed in the HI horizon (0—0.5 m) 
enhanced precise data collection, resulting in better measurement accuracy than in the H, horizon (0.5-2.0 m). The study also 
discussed the factors responsible for the deviation between measured SWC and estimated SWC. The development and 
adoption of precision farming and rational irrigation required detail knowledge of soil and crop. The method used in this study 
was useful in the research and description of high-resolution soil spatial variability and hydric characteristics. Additional field 
calibration was required for applying the method practical on routine field application. The use of general petro-physical 


relationship between soil electrical resistivity and its moisture, if appropriate, could bring the method a step closer to practical 


field application for the purpose of improving irrigation management. 
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在 干旱 半 干 旱 区 土壤 含水 量 是 影响 农作物 生长 
的 重要 因素 ， 控 制 和 监测 土壤 水 分 变化 规律 ， 对 农 
业 生态 环境 有 着 重要 的 意义 10。 过 去 几 十 年 中 ,很 
多 研究 人 员 已 经 尝试 多 种 方法 来 监测 土壤 合 水 量变 
化 规律 。 中 子 水 分 仪 是 一 种 测量 深层 土壤 含水 量 与 
水 分 动态 变化 的 有 效 方 法 , 但 其 放射 源 受 众多 约束 
条 件 限 制 咎 ,时 域 反 射 仪 和 电容 探 针 是 近年 来 兴起 
的 一 种 实时 测量 浅 层 土 壤 合 水 量 的 方法 ?这些 
方法 仅仅 能 提供 距离 传感器 探头 较 近 区 域 的 土壤 合 
水 量 信息 ， 大 尺度 范围 应 用 效果 较 差 ， 费 用 高 ; TH 
反 卫 星 和 遥感 探测 技术 能 大 范围 监测 土壤 水 的 空间 
分 布 , 但 是 探测 深度 受 限于 几 厘 米 ， 空 间 分 布 率 较 
低 , 并 且 受 植被 、 微 地 貌 等 因素 的 影响 中。 

相对 而 言 ， 电 磁 、 探 测 雷 达 、 电 阻 率 等 地 球 物 
理 方法 测量 土壤 含水 量 不 影响 土壤 结构 ， 测 量 结 
能 覆盖 地 球 表层 的 空间 变异 性 , 测量 范围 能 达 数 十 
米 其 至 上 百 米 。 其 中 ,电阻 率 成 像 技术 已 经 被 广泛 
运用 到 水 文 地 质 探 测 中 ”近年 来 ,这 种 方法 也 
用 来 研究 土壤 表层 降水 入 渗 与 土壤 含水 量 25 4、 潜 
在 补给 5、 植物 根系 时 空 分 布 规律 几 和 地 下 水 动 
SM, Michot 等 名 研究 表明 高 密度 电阻 率 成 像 法 
(ERT) 能 界定 土屋， 以 便于 监控 灌溉 期 玉米 (Zea 
mays) 土 壤 水 运 移 ; Wenninger 等 "和 Koch 等 ("通过 
结合 水 文 经 典 测 量 和 示 踪 技术 ,研究 发 现 ERT 可 以 
识别 坡 面 水 流 路 径 ; Schwartz 等 用 二 维 电 阻 率 图 
像 、 一 维 TDR 土壤 合 水 量 测量 结合 土壤 物理 化 学 特 
性 建立 了 非 均 质 黏 性 土壤 的 二 维 含水 量 模型 。 这 些 
都 为 建立 干旱 半 干 旱 区 土壤 含水 量 与 电阻 率 相关 关 
系 模型 贷 定 了 基础 。 然 而 目前 国内 应 用 ERT 对 干旱 


Electrical resistivity tomography; Soil water content; Rainfall infiltration; Unsaturated soil; Time domain 


半 干 旱地 区 土壤 合 水 量 监测 研究 成 果 很 少 。 本 文 以 
陇 中 半 干 旱 区 玉米 耕作 地 为 例 ， 通 过 在 土壤 表面 布 
EEI, 分 别 在 降雨 前 和 降雨 后 进行 二 维 ERT 数据 
的 测量 ， 实 现 对 土壤 二 维 剖 面 含 水 量 和 电阻 率 值 监 
测 , 分 析 不 同 条 件 下 土壤 含水 量变 化 的 原因 ， 建 立 
陇 中 半 干 旱 区 农田 土壤 电阻 率 和 合 水 量 之 间 的 相关 
关系 ,并 用 实测 的 ERT 数据 估计 土壤 合 水 量 , 为 实 
现 精 确 和 高 效 的 农业 管理 提供 一 种 新 途径 。 


1 研究 区 概况 与 研究 方法 
11 ”研究 区 概况 

试验 田 位 于 兰州 大 学 半 干 旱 区 农业 生态 系统 试 
验 站 (35.42°N, 105.24°E), 海拔 高 度 1951 m 研究 区 
气候 属于 典型 的 温暖 带 半 干 旱季 风气 候 ,， 年 平均 降 
水 量 350 mm， 主 要 集中 在 7、8、9 月 ,年 平均 蒸发 
量 1 800 mm, 无 霜 期 约 136~186 d， 屋 夜 温 差 大 ， 光 
照 充足 。 试验 田 为 常规 种 植 农田 ， 面 积 为 20 mx30 m, 
试验 期 间 种 植 玉米 ， 玉 米 行距 和 株距 分 别 为 0.5 m 
和 0.4 mo 
1.2 ”试验 场地 

为 了 能 够 有 效 地 采集 试验 田 二 维 剖 面 ERT 成 像 
数据 , 按照 试验 田 纵 长 方向 , 以 0.5 m 为 间距 (如 图 1 
所 示 ), 布置 了 一 列 52 个 电极 ,电极 埋 深 为 0.5 mo 
每 个 电极 依次 与 多 芯 电缆 自动 转换 开关 一 一 相 接 。 
13 ”试验 准备 阶段 

不 同 深度 土壤 含水 量 和 温度 数据 用 Stevens 
Water Monitoring Systems 采集 。 其 中 土壤 含水 量 的 
监测 主要 用 HydraProbe I (A FEE HP 1) 探 头 ， 
基于 时 域 反 射 原理 ， 为 了 便于 监测 土壤 含水 量变 化 
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(a 和 为 实际 电极 排列 , e 为 校准 坑 Q BUE, d 为 电极 、 含 水 量 探头 分 布 ) 


Fig. 1 


Spatial distribution maps of electrodes and soil water content monitoring in the field site for imaging data of high-density 


electrical resistivity tomography (ERT) (a, b show real electrodes arrangement; c is soil section of calibration pit Q; d shows the 


distribution of electrodes and HydraProbe ll ) 


和 建立 含水 量 与 电阻 率 之 间 的 关系 , 在 试验 场地 内 ， 
控 4 DEH M, N, P, Q, 如 图 1 所 示 。 其 中 基 坑 
M、N、P 用 来 实时 监测 土壤 水 分 和 温度 变化 , 在 2m 
的 深度 内 ,将 基 坑 M、N 和 P 的 8 个 垂 向 深度 (0.1 m. 
0.2m, 0.3m, 0.5m, 0.8m, 1.0m 1.5m, 2.0m) 
分 别 安装 相同 的 HP 1 探头 ， 然 后 再 将 基 坑 一 层 一 
层 填 埋 , 为 了 防止 土壤 结构 破坏 对 ERT 数据 产生 
影响 ,对 3 个 基 坑 提前 6 个 月 进行 HP1I 的 安装 。 基 
H Q 是 在 ERT 数据 测量 时 段 前 1 d 进行 挖掘 ， 
Michot 等 ?HJ 证 实 对 土壤 进行 分 层 考虑 可 以 更 为 准 
确 地 建立 土壤 含水 量 与 电阻 率 值 之 间 的 关系 ,根据 
土壤 颗粒 物理 和 化 学 性 质 ， 把 土壤 分 为 两 层 0~0.5 m 
(HDA 0.5-2.0 m (H2) 分 别 进行 土壤 电阻 率 的 测定 。 
BOEN Q 电阻 率 测量 采用 Wenner 的 排列 方式 ， 电 
极 间距 为 0.2 m， 埋 深 为 0.5 m， 水 平 排列 在 土 层 剖 
E 0.25 m 和 1.25 m 的 深度 处 (如 图 2 所 示 ); 另外 
校准 坑 Q 的 土壤 含水 量 和 温度 分 别 用 HP 1 探头 和 
铂金 丝 温度 探 针 测 定 ， 对 HPI 探头 的 矫正 采用 
Heimovaarat 提出 的 方法 。 表 1 为 H 和 H, AE 
塘 理 化 性 质 。 

1.4 ERT 数据 采集 过 程 


本 次 试验 采集 仪器 使 用 西安 澳 立 华 勘 探 技术 开 


发 有 限 公司 生产 的 FlashRES-UNIVERSAL 多 通道 、 
超 高 密度 直流 电 法 勘探 系统 ，ERT 数据 的 采集 也 采 
用 Wenner 排列 方式 。 如 图 3 所 示 ， 由 数据 采集 发 射 
机 产生 一 个 恒定 的 输入 电流 , 通过 指定 的 电流 电极 
对 (CiC?) 输 入 到 测量 介质 内 部 ,形成 一 个 人 工 电流 
场 。 同 时 在 计算 机 的 控制 下 ,接收 机 根据 事先 给 定 
的 测量 组 合 指定 的 测量 电位 电极 对 (P1P;)， 测 量 相 
应 的 电位 差 以 及 电流 的 大 小 ， 并 保存 测量 结果 ， 如 
此 循环 直到 测量 完 所 有 事先 给 定 的 电极 测量 组 合 
为 止 。 
测量 的 电位 或 电位 差 从 一 定 程度 上 反映 了 测量 
介质 电导 率 的 变化 信息 。 利 用 测量 的 电位 或 电位 差 ， 
通过 成 像 算法 可 重建 出 场 内 介质 的 电导 率 的 空间 分 
布 特征 ， 进 而 进行 物性 特征 的 分 析 与 解释 ， 实现 了 
可 视 化 成 像 测量 。Wenner 采集 方式 对 于 垂 向 电阻 率 
变化 非常 敏感 ， 有 利于 得 到 更 准确 的 反 演 结果 。 介 
质 视 电 阻 率 o 值 计 算 采 用 如 下 公式 : 
p 2 2na(AV / T) (1) 
式 中 : 7 为 电极 Cl 和 C; 之 间 的 输入 电流 , AV 73 Pi、 
P5 两 个 电极 之 间 的 电位 差 , a 为 电极 间距 。 
试验 的 测量 时 段 是 2015 年 8 月 20—26 H, 
进行 8 次 ERT 数据 的 采集 ,如 图 4 所 示 。 其 中 8 月 


+ 
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电极 
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2 校准 坑 Q 不 同 层次 土壤 地 球 物理 特性 建立 示意 图 
Fig.2 Experimental setup for geophysical characterization of soil horizons in the calibration pit Q 
表 1 试验 场地 不 同 层次 土壤 颗粒 物理 和 化 学 性 质 
Table 1 Physical and chemical properties of soil particle in different horizons of field site 
dr y My Dax 7S 
n 埋 深 ioi 粉 粒 Esp caco， 干 容重 
dodien Depth (m) (0.06 92) (0.002 x d«-0.06) (d«0.002) (96) Bulk density p 
Sand (%) Silt (96) Clay (96) (gem?) 
Hi 0~0.5 27,5 47.9 9.8 4.0 1.35 8.1 
H2 0.5-2.0 38.3 31.6 15.1 25.8 1.52 8.3 


Dd 为 土壤 颗粒 粒 径 ， 单 位 为 mm。1) d is diameter of soil particle, whose unit is mm. 
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; 位 置 标识 点 
n=3 十 十 十 十 十 二 十 十 Signature point of 
CiCs 电 流 电极 PP, 电 压 电 极 data space position 


Current electrode Voltage electrode 


图 3 高 密度 电阻 率 成 像 法 (ERT) 数 据 采 集 方法 (Wenner 排列 ) 和 电阻 率 空间 分 布 图 


Fig.3 Schematic map of high-density electrical resistivity tomography (ERT) data measurement (using Wenner electrode array) and 


the spatial distribution map of electrical resistivity 


22 日 17:00 试验 田 下 了 一 场 中 雨 ,持续 了 大 约 3.5 ho 实 的 土壤 电阻 值 ， 所 以 进行 ERT 数据 反 演 显得 尤为 


前 3 次 ERT 数据 采集 是 在 降雨 前 ,其 余 5 次 均 在 降 重要。 电阻 率 的 反 演 计算 属于 位 场 反 演 王 )， 


即 通过 


雨 结束 之 后 。 下 十 期间 未 进行 ERT 数据 采集 ， 所 以 一 定 的 反 演 方法 寻找 出 一 个 能 够 给 出 与 测量 值 具有 


雨水 对 ERT 数据 准确 性 影响 很 小 。 相同 响应 特征 的 参数 化 模型 ， 即 能 用 有 限 参 数 表征 
15 ERT 数据 反 演 的 模型 ， 该 模型 是 测量 介质 在 数学 上 理想 化 的 一 个 


由 公式 (1) 实 测 的 视 电 阻 率 数 据 本 身 不 能 给 出 真 ” 再 现 "1 在 进行 数据 反 演 之 前 , ERT 采集 的 数据 需 


http://www.ecoagri.ac.cn 


第 10 期 后“ 宁 等 ; 基于 高 密度 电阻 率 成 像 法 的 陇 中 半 干 时 区 土壤 含水 量 监测 研究 in 了 


-H1421 


o ERT 数据 测量 编号 
Measuring number 
dd of ERT data 


降水 量 
Precipitation (mm) 
D 


Z 合 
IE ARRE 


1 1 | | 
Eee bnc Le eed 


08-22 08-23 08-24 08-25 08-26 
测量 日 期 (月 -日 ) 和 时 间 


Date (month-day) andtime ofmeasurement 


4 高 密度 电阻 率 成 像 法 (ERT) 数 据 采 集 和 自然 降水 时 
间 记 录 
Fig.4 Times of high-density electrical resistivity tomography 
(ERT) data acquisition and natural rainfall records 


检查 ， 并 消除 受 噪声 影响 较 大 的 数据 。 电 流 输 入 小 
于 0.01 A 或 者 是 相对 误差 大 于 5% 的 数据 在 反 演 过 
程 中 应 删除 ?23H。Panissod 等 1 表明 二 维 电 阻 率 反 
演 相 上 比较 三 维 反 演 在 反映 剖面 土壤 合 水 量 时 更 为 恰 
当 。 所 以 , 为 了 获得 真实 土壤 剖面 电阻 率 分 布 ， 本 文 
采用 RES2DINYV 软件 反 演 计算 ERT 数据 中 ,该 软件 
基于 圆滑 约束 条 件 下 的 最 小 二 乘 方法 ， 建 立 二 维 的 
土壤 电阻 率 剖 面 模型 ， 圆 滑 约 束 最 小 二 乘法 基于 以 
下 方程 : 


08-20 08-21 


(JJ *uF)d 2 J'g (2) 

式 中 : FEAR 了, 大 为 水 平平 滑 滤 波 系 数 和 矩阵 , f. 73 
垂直 平滑 滤波 系数 和 矩阵 , 7 为 偏 导 数 和 矩阵 , .7 为 /的 转 
iBABEE, u 为 阻尼 系数 , 4 为 模型 参数 修改 矢量 , g 为 
残 差 矢量 。 

RES2DINV 软件 第 1 次 迭代 过 程 中 ,用 电阻 率 
偏 导数 计算 的 解析 值 建立 的 均匀 土壤 电阻 模型 当 作 
初始 模型 。 随 后 的 迭代 过 程 ， 运 用 拟 牛 顿 法 估计 偏 
导数 , 减少 了 电脑 计算 时 间 。 反 演 模型 由 和 矩形 网 格 
构成 ,通过 对 区 域 ERT 数据 测量 ,， RES2DINV 软件 
能 给 出 一 个 反 演 计算 的 二 维 电 阻 率 剖 面 ， 迭 代 优 化 
的 方法 能 减少 实测 电阻 率 值 和 模型 反 演 计算 得 到 的 
电阻 率 值 之 间 的 误差 值 (RMS, %)。 
1.6 ”土壤 合 水 量 电 阻 率 特性 的 确定 

土壤 电阻 率 大 小 往往 受 土壤 结构 、 孔 隙 水 电阻 、 
土壤 舍 水 量 、 土 壤 温 度 等 一 系列 因素 影响 3。 测量 
期 间 ， 雨 水 和 地 下 水 的 电导 率 展现 出 很 微小 的 变化 
[(0.033+0.007) S.m 和 (0.029+0.005) S.m ']， 表 明 孔 
隙 水 的 电导 率 相对 于 土壤 合 水 量 的 变化 较为 恒定 。 
因此 本 文 和 其 他 ERT 研究 一 样 1， 可 以 在 土壤 电 
阻 率 测量 中 忽略 孔隙 水 电阻 率 的 影响 。 

电阻 率 在 不 同 温 度 和 土壤 含水 量 条 件 下 的 实测 
值 不 同 。 为 了 消除 温度 的 影响 ， 确 定 土壤 电阻 率 和 


含水 量 之 间 的 单一 关系 , 所 有 的 测量 结果 都 要 在 
25 标准 温度 下 矫正 。 本 文 假定 试验 过 程 中 水 平 剖 
面 温度 分 布 均 匀 , 在 0~1.0 m 深度 上 满足 线性 分 布 ， 
温度 矫正 采用 25 CH FHI Keller and Frischknecht 
AATRES, 
Pref = Psoil [1+a(Ton -T«)] (3) 

式 中 : ou (Qm) 是 在 25 CERE TA('C) 条 件 下 
的 电阻 矫正 值 ， Psoii (Qm) 是 在 土壤 温度 Toil CR 
下 反 演 的 电阻 率 值 ; o 是 修正 系数 ,一 般 等 于 0.025。 

研究 通过 线性 回归 方法 采用 校准 坑 Q 实测 的 电 
阻 值 和 土壤 含水 量 数据 来 建立 不 同 土壤 分 层 (H, 和 
H2) 条 件 下 两 者 之 间 的 相关 关系 ， 如 图 5 所 示 。 
Michot 等 J 研究 证 实 多 项 式 画 数 或 者 害 画 数 必须 运 
用 于 完全 饱和 状态 和 干 土 状 态 之 间 的 大 范围 的 土壤 
水 分 变化 ,然而 对 于 自然 条 件 下 体积 含水 量 在 永久 
性 鞭 戎 点 与 田间 持 水 量 之 间 的 情况 下 ， 约 10% 到 
25%, 一 元 线性 函 数 适 用 性 很 好 。 
2 利用 ERT 数据 计算 土壤 含水 量 

本 文通 过 二 维 ERT 数据 建立 矩形 网 格 来 反 演 计 
算得 到 二 维 土壤 含水 量 剖 面 , 如 图 5 所 示 ， 和 矩形 网 
格 由 385 个 小 的 矩形 单元 组 成 。 中 心 点 的 坐标 定义 
73 (X, Z), X RE BAR BERI 75 TBSZKSE RB S, Z RRE 
形 模型 块 中 心 点 的 垂 向 深度 。 网 格 剖 分 为 12 2, 各 
层 深 度 分 别 为 0.25 m、0.50 m、0.775 m、1.08 m、 
1.41 m. 1.78 m, 2.18 m, 2.62 m, 3.11 m, 3.64 m 
和 4.23 mo 

对 于 每 个 矩形 模型 单元 ， 中 心 土壤 温度 用 线性 
差分 计算 ， 中 心 点 土壤 含水 量 通过 中 心 点 的 土壤 电 
阻 率 值 计算 得 到 ,三 者 之 间 矫 正 关 系 符合 垂直 方向 
真实 的 土壤 温度 梯度 。 二 维 土壤 含水 量 剖 面 通过 模 
型 单元 中 心 点 的 含水 量 值 三 角形 剖 分 得 到 。 为 了 验 
证 土壤 合 水 量 预测 的 准确 性 , 用 SWMS RE 4 次 
(R2、R4、R6、R8)M、N、P 基 坑 不 同位 置 处 的 土壤 
合 水 量 值 。 
3 结果 和 分 析 
31 土壤 含水 量 和 温度 分 析 

图 6 给 出 了 一 次 降水 过 程 4 次 测量 时 间 (R2 为 
降水 前 , R4、R6 和 R8 为 降水 后 ) 实 测 土壤 体积 含水 
量 和 校准 坑 Q 温度 变化 情况 。 从 图 中 可 以 出 , 由 于 
降水 的 影响 ,在 0~1.0 m 深度 范围 内 ,土壤 含水 量 都 
有 了 明显 增加 ， 位 置 不 同 , 含水 量变 化 情况 也 不 同 。 
N、P 剖面 处 土壤 透水 性 较 好 , 0~1.0 m REALES 
水 量 明显 增加 ， 相 比较 M 剖面 处 含水 量变 化 不 大 ， 
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5 不 同 分 层 条 件 下 土壤 剖面 含水 量 和 电阻 率 之 间 的 矫正 关系 ( 皇 2$S 'C) 以 及 利用 电阻 率 计 算 二 维 体积 含水 量 剖 面 过 
程 示意 图 


Fig. 5 Calibration (1-25 'C) between soil volumetric water content and electrical resistivity of different horizons and the schematic 
diagram of calculating 2D volumetric water content of soil by using electrical resistivity 
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6 4 次 高 密度 电阻 率 成 像 法 (ERT) 采 集 时 间 [R2( 降 雨 前 ) 和 R4、R6、R8( 降 水 后 )] 的 剖面 M、N 和 了 实测 土壤 体积 含 
水 量 和 Q 剖面 实测 土壤 温度 变化 
Fig.6 Variations of measured volumetric water contents of M, N, P profiles and temperature of Q profile during four various 
high-density electrical resistivity tomography (ERT) surveys [R2 (before rain) and R4, R6, R8 (after rain)] 
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可 能 与 该 处 土壤 压 实 情况 有 关 。 降 水 之 前 , 土壤 含水 
量 最 大 值 主要 集中 在 P 坑 深度 为 0.2 mAb, 最 小 值 集中 
在 P 坑 深 度 为 1.0 m 处 ,实测 土壤 含水 量 最 大 值 为 
27.6%， 最 小 值 为 12.4%。 降 水 之 后 , 0-0.5 m 深度 处 土 
壤 含 水 量 增长 幅度 最 大 , 土壤 含水 量 最 大 值 主要 集中 
在 P 坑 深度 为 0.2 m 处 , 最 小 值 集中 在 P 坑 深 度 为 0.8 m 
处 , 实测 土壤 含水 量 最 大 值 为 30.2%, 最 小 值 为 12.1%。 

土壤 温度 变化 和 土壤 含水 量变 化 之 间 并 没有 直 
接 的 相关 关系 。 降 水 之 前 ,土壤 温度 随 着 埋 深 的 增 
加 大 体 上 呈 下 降 的 趋势 ， 降 水 之 后 ,在 0~0.5 m 深度 
内 由 于 昼夜 变化 和 降水 等 因素 的 影响 ,土壤 表层 温 
度 下 降 明 显 ， 而 在 0.5~1.0 m 深度 处 ,温度 有 了 略微 的 
回升 ， 满 足 前 文 线 性 分 布 的 假定 , 1.0 m 深 度 以 下 温度 
基本 趋 于 稳定 。 因此 在 ERT 数据 获取 的 过 程 中 , 日 周 
期 和 水 分 入 渗 是 影响 土壤 温度 变化 的 重要 因素 1。 
3.22 ”基于 二 维 土壤 剖面 电阻 率 图 像 分 析 

在 ERT 处 理 过 程 中 ,， 本文 对 4 次 ERT(R2、R4、 
R6、R8) 测 量 数据 进行 反 演 计算 ， 得 到 4 次 测量 的 土 
塘 剖 面 二 维 电阻 率 图 像 (图 7), 4 次 测量 数据 生成 的 
图 像 分 别 用 GI. G2. G3. G4 m. Él G1 是 试验 
开始 测量 时 土壤 电阻 率 反 演 结果 ， 在 降水 入 渗 之 前 ， 
表明 初始 的 土壤 结构 和 电阻 率 值 分 布 特征 ; 图 G2、 
G3. G4 是 降水 之 后 测量 土壤 电阻 率 值 的 反 演 结果 ， 
表明 随 着 降水 入 渗 过 程 的 进行 , 土壤 电阻 率 值 的 变 
化 情况 。 从 图 中 可 以 看 出 : 

1) 不 同时 刻 测量 数据 获得 的 土壤 电阻 率 图 像 在 
二 维 平面 中 的 分 布 是 不 同 的 ， 相 对 于 图 G1 fmm, 
水 入 渗 使 得 二 维 剖 面 土壤 电阻 值 整体 上 呈现 出 明显 
降低 的 趋势 ,尤其 是 在 土壤 垂 向 深度 为 0.2-1.0 m 和 
3.0-4.5 m 范围 内 的 电阻 率 值 相对 较 低 。 

2) 对 于 垂 向 深度 为 0~0.2 m 的 区 域 , 土壤 电阻 率 
值 相对 较 高 ， 一 方面 由 于 试验 场地 常年 种 植 玉米 ， 
土壤 表层 覆盖 有 地 膜 ， 对 表层 土壤 有 保护 作用 ; 5 
一 方面 由 于 玉米 根系 的 固 结 作用 和 土壤 压 实 的 影响 ， 
Ankeny 等 中 "证 实 了 压 实 作 用 所 造成 的 土壤 大 孔 阶 
的 破坏 会 降低 土壤 的 渗透 系数 ， 致 使 降水 不 能 均匀 
入 渗 到 整个 试验 的 平面 内 。 但 在 一 些 倾 斜 的 坡 面 处 
或 者 是 坑 洼 处 ， 降 水 汇集 较 多 ， 下 渗 明 显 ， 即 选择 
性 入 渗 现 象 。 例 如 图 G2 Fh, 在 X-4.0-5.0 m, RE 
为 0~0.2 m, 存在 着 电阻 率 较 低 的 异常 值 。 

3) 在 0.2~1.0 m 的 垂 向 深度 内 , G2、G3 存在 着 很 
多 不 符合 降水 入 渗 局 部 异常 区 域 ， 例 如 图 G2 和 G3 
在 X-9.0-10.0 m 处 。 本 应 该 电阻 率 值 降低 的 地 方 ， 
却 反而 增 大 。 这 种 异常 通常 是 由 于 玉米 根部 对 局 部 
土壤 水 分 的 吸收 、 土 壤 质 地 变异 性 和 入 渗 的 非 均 匀 


性 (优先 流 的 存在 ) 等 因素 造成 的 I， 导致 局 部 的 土 
塘 含 水 量 反而 降低 ， 土 壤 电阻 率 值 增 大 。 即 在 土壤 
湿润 阶段 ， 同 一 位 置 不 同时 间 点 的 土壤 电阻 率 值 变 
化 规律 也 不 相同 ， 存 在 着 阻 值 异 常 增 大 的 情况 。 

4)G1 图 中 在 X-7.0-11.0 m 和 X-13.0-19.0 m, 
深度 为 1.5-3.0 m 的 范围 内 ,存在 两 块 明 显 的 高 阻 值 
区 域 ， 由 于 降水 下 渗 作 用 的 影响 G2、G3 图 中 该 区 
域 电 阻 率 值 明显 降低 ，G4 图 中 电阻 率 值 反而 增高 ， 
这 种 电阻 率 值 “高 - 低 -高 ”的 变化 过 程 , 与 一 次 降水 
过 程 前 后 “ 干 - 湿 - 干 ” 的 循环 过 程 颇 为 一 致 ， 体 现 出 
ERT 图 像 能 定性 地 反映 出 含水 量变 化 趋势 。 
3.3 ERT 估计 土壤 含水 量 与 实测 土壤 含水 量 比较 

根据 二 维 剖 面 土壤 电阻 值 、 土 壤 仿 水量、 土壤 
温度 三 者 之 间 的 矫正 关系 ,对 4 次 测量 (R2、R4、R6、 
R8) 的 ERT 数据 进行 处 理 ， 计 算 二 维 剖面 的 土壤 合 
水 量 数据 ， 反 演 得 到 4 个 二 维 剖 面 土壤 含水 量 图 像 ， 
对 应 为 W1、W2、W3、W4( 图 8)。 从 图 中 可 以 看 出 ， 
二 维 剖 面 土壤 含水 量 的 空间 结构 与 电阻 率 剖 面 的 空 
间 结 构 是 相似 的 ， 垂 向 入 渗水 流 路 径 的 变化 符合 电 
阻 率 的 变化 , 并 且 可 以 明显 地 看 出 降水 入 渗 “ 干 - 湿 - 
干 ” 循 环 过 程 。 

为 了 证 明 建 立 的 ERT 估计 土壤 含水 量 的 准确 性 ， 
本 文 利用 (M、N、P) 3 个 基 坑 实测 的 土壤 含水 量 数 
据 和 ERT 估计 的 土壤 含水 量 数据 建立 了 两 者 之 间 的 
相关 关系 (如 图 9， 表 2)。 可 以 看 出 ,实测 的 土壤 含水 
量 值 与 估计 的 土壤 含水 量 值 之 间 有 一 个 较为 显著 的 
线性 关系 (R=0.651 8,，n=96)， 实测 值 与 估计 值 的 余 
率 近 似 等 于 1， 用 均 方 根 误差 (RMSE, %) 代 表 估 计 精 
度 ， 用 平均 误差 (ME, %) 代 表 估 计 人 和 偏差。 通过 表 2 看 
tH, 在 0~2.0 m 深度 总 体 上 估计 偏差 较 小 ， 为 0.74%; 
土壤 含水 率 的 估计 精度 较 高 ， 为 2.67%。 在 0~0.5 m 
深度 范围 内 ，HP II 探头 分 布 较为 密集 ， 数 据 采 集 较 
为 准确 , 故 Hi 层 的 估计 精度 略 高 于 Ho 层 。 本 文 提 供 
了 一 个 高 分 辩 率 的 土壤 结构 二 维 分 布 与 水 分 运 移 过 
程 的 图 像 ， 利 用 土壤 电阻 率 和 含水 量 之 间 建 立 的 岩 
石 物 理 关 系 ， 可 以 定性 的 反映 土壤 二 维 剖 面 土 壤 合 
水 量 的 变化 过 程 ， 但 是 介 于 电极 间距 布置 较 大 ， 很 
难 精确 的 估计 出 单个 位 置 处 土壤 含水 量变 化 。 

本 文通 过 野外 试验 发 现 ， 影 响 实测 的 土壤 含水 
量 和 ERT 估计 的 土壤 含水 量 之 间 偏 差 的 因素 有 很 多 ， 
主要 体现 在 : 1) 相 比较 小 尺度 的 土壤 水 分 变化 过 程 ， 
本 试验 安装 的 HP 1 探头 之 间 的 间距 要 大 得 多 , 不 能 
很 好 地 反映 土壤 含水 量 的 变化 过 程 。2) 实 测 的 土壤 
电阻 率 数据 和 本 文 反 演 过 程 得 到 的 电阻 率 数 据 不 同 ， 
电极 的 布置 、 反 演 精 度 、 电 极 间距 等 一 系列 不 确定 
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二 维 土 壤 剖 面 电阻 率 随 时 间 变 化 过 程 [4 次 高 密度 电阻 率 成 像 法 (ERT) 采 集 , G1( 降 水 前 ) 和 G2、G3、G4( 降 水 后 )] 


Fig.7  Two-dimension maps of electrical resistivity temporal variations of soil profiles during four various high-density electrical 


图 8 
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resistivity tomography (ERT) surveys [G1 (before rain), G2, G3, and G4 (after rain)] 
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二 维 土壤 剖面 体积 含水 量 随时 间 变 化 过 程 [4 次 高 密度 电阻 率 成 像 法 (ERT) 采 集 ，W1( 降 水 前 ) 和 W2、W3、W4( 降 水 后 )] 


Fig.8 Two-dimension maps of soil volumetric water content temporal variations of oil profiles during four various high-density 


electrical resistivity tomography (ERT) surveys [W1 (before rain), W2, W3, and W4 (after rain)] 


http://www.ecoagri.ac.cn 


第 10 期 E TE 基于 高 密度 电阻 率 成 像 法 的 陇 中 半 干 旱 区 土壤 含水 量 监测 研究 上 


-H1425 


30 F 深度 Depth (m): k P 
274 0.1; 0.2; 0.3; 0.5 A 
LT. 
24L »720.9397x12.357 A A 
R^-0.7350 
21 
18 F 
15 F 
一 一 一 95% 置 信 水 平 
12 要 Confidence level 
o 1 1 1 1 1 1 
30r 深度 ; ^ 
深度 Depth (m): 
0.8; 1.0; 1.5; 2.0 í pg 


25. y=0.940 3x-0.174 
R^-0.683 8 AA, 
N-48 全 
20 


实测 土壤 含水 量 
Measured moisture (96) 


35r 深度 Depth (m): 
0.1; 0.2; 0.3; 0.5 
30L 0.8;1.0; 1.5; 2.0 A 


y70.911 0x+1.735 
R^-0.6518 


L 
9 12 15 18 21 24 2T 30 
ERT 估计 土壤 含水 量 Estimated moisture(96) 


9 不 同 层次 高 密度 电阻 率 成 像 法 ERT) 估 计 的 土壤 含 
水 量 与 实测 土壤 含水 量 相关 关系 分 析 
Fig.9 Correlation analysis between estimated soil moisture by 


high-density electrical resistivity tomography (ERT) and 
measured soil moisture of different horizons 


的 因素 都 会 导致 两 者 之 间 的 偏差 。3) 尽 管 小 电极 可 
以 减少 电极 和 土壤 之 间 的 接触 电阻 , 但 是 对 于 浅 层 


探查 深度 和 高 分 辨 率 的 试验 要 求 来 说 , 把 电极 考虑 
成 点 源 来 处 理 , 还 是 非常 困难 的 。 4) 土 壤 电 阻 率 变化 
一 般 会 受到 降水 入 渗 灌 后 现象 的 影响 。5) 土 壤 的 传 
导电 流 主要 是 由 电极 发 出 的 , 影响 土壤 电阻 率 值 的 
大 小 除了 土壤 湿度 和 温度 之 外 ,还 会 与 土壤 中 各 种 
离子 的 分 布 浓度 有 关 。 降 水 的 不 均匀 入 渗 会 导致 局 
部 地 区 土壤 溶液 的 TDS 发 生变 化 ,造成 了 不 可 忽视 
的 土壤 电阻 率 的 变化 :。 以 上 这 些 方面 的 问题 ， 还 
需要 以 后 继续 研究 。 
4 结论 

近年 来 ,用 地 球 物 理 测量 的 方法 研究 高 分 辩 率 
的 地 下 水 流 的 流动 过 程 ， 由 于 其 本 身 不 损害 土壤 微 
结构 ,而 且 测 量 方法 不 繁 珊 ， 已 经 变 得 越 来 越 受 到 
欢迎 和 重视 。 本文 以 陇 中 半 干 旱 区 玉米 耕作 地 为 例 ， 
通过 在 土壤 表面 布置 电极 , 用 二 维 ERT 数据 和 土 
塘 仿 水量 数据 ,建立 了 陇 中 半 干 旱 区 土壤 电阻 率 
和 合 水 量 之 间 的 相关 关系 ,通过 对 一 场 降水 前 后 
土壤 含水 量 的 变化 ,实现 对 二 维 剖面 土壤 含水 量 
的 动态 监测 。 试 验 结果 表明 , 土壤 合 水 量 的 变化 可 
以 通过 土壤 电阻 率 变化 体现 出 来 ， 二 维 的 土壤 合 
水 量 剖 面 通过 二 维 ERT 数据 和 两 者 之 间 建 立 的 相 
关 关 系 反 演 计算 得 到 ,土壤 温度 的 变化 在 计算 中 同 
样 被 考虑 到 。 

土壤 合 水 量 在 一 次 降水 前 后 有 较为 明显 的 “ 干 - 
湿 - 干 ”循环 变化 过 程 。 降 水 之 前 , 气候 较为 干燥 , E 
壤 电 阻 率 整 体 偏 高 ,降水 入 渗 使 得 二 维 剖 面 土壤 电 
阻 率 整体 上 呈现 出 明显 降低 的 趋势 。 由 于 土壤 压 实 
作用 的 影响 ,表层 土壤 电阻 率 变 化 不 明显 , 但 是 在 
某 些 土质 相对 松散 的 局 部 区 域 ， 也 存在 着 合 水 量 的 
大 量 富 集 ， 入 渗 过 程 中 , 受 土壤 质地 变异 性 、 入 渗 的 
非 均 匀 性 (优势 流 )、 玉 米 根部 对 土壤 水 分 的 吸收 等 因 
素 影响 ， 也 会 导致 局 部 区 域 含水 量 反 而 更 低 。 


R2 高 密度 电阻 率 成 像 法 (ERT) 估 计 土 壤 含 水 量 与 实测 土壤 含水 量 统计 分 析 结 果 
Table 2 Statistical analysis results between estimated soil moisture by high-density electrical resistivity tomography (ERT) and 
measured soil moisture 


3878 Depth (m) n R R? ME (%) RMSE (%) 
Hı: 0.1;0.2;0.3;0.5 48 0.8573 0.735 0 1.28 2.58 
H»: 0.8;1.0;1.5;2.0 48 0.826 9 0.683 8 -2.02 2.70 

0.1 ~ 2.0 96 0.807 3 0.6518 -0.74 2.64 


实测 的 土壤 含水 量 值 与 估计 的 土壤 含水 量 值 之 
间 有 一 个 较为 显著 的 线性 关系 (R =0.651 8, n-96), 
实测 值 与 估计 值 的 斜率 近似 等 于 1。 在 0-2.0 m 深度 
范围 内 ,总 体 上 估计 偏差 较 小 , 为 0.74%; 土壤 合 水 


率 的 估计 精度 较 高 , 为 2.64%， 在 0~0.5m 深度 范围 
内 ，HP 1 探头 分 布 较为 密集 ， 数 据 采 集 较为 准确 ， 
故 H 层 的 估计 精度 略 高 于 H 层 。 但 是 本 文 电极 间 
距 布 置 较 大 ,很 难 精确 地 估计 出 单个 位 置 处 土壤 合 
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水 量变 化 过 程 。 [10] Binley A, Kemna A. DC resistivity and induced polarization 
当今 精 准 农业 和 合理 灌溉 需 要 更 为 成 熟 的 作物 methods[M]//Rubin Y, Hubbard S S. Hydrogeophysics. 
和 土壤 知识 体系 , 相 比 较 之 前 利用 实测 工具 进行 野 AE a penu NIE 
à : R Riu [1] 苏 永 军 , 王 绪 本 ， 罗 建 群 ， 高 密度 电阻 率 法 在 三 星 堆 壕沟 
外 测量 ， 地 球 物理 方法 精度 较 高 。 本 文 提 供 了 一 考古 勘探 中 应 用 研究 []. 地 球 物理 学 进展 ，2007,，22(1): 
高 分 辩 率 的 土壤 结构 二 维 分 布 与 水 分 运 移 过 程 的 图 268-272 
像 ， 利用 土壤 电阻 率 和 水 分 之 间 建 立 的 岩石 物理 关 Su Y J, Wang X B, Luo J Q. The ahaeological application of 
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: " Sanxingdui Site[J]. Progress in Geophysics, 2007, 22(1): 
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